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Abstract. We consider the problem to reconstruct a wave speed 
c e C°°(M) in a domain M C W l from acoustic boundary mea- 
surements modelled by the hyperbolic Dirichlct-to-Neumann map 
A. We introduce a reconstruction formula for c that is based on 
the Boundary Control method and incorporates features also from 
the complex geometric optics solutions approach. Moreover, we 
show that the reconstruction formula is locally Lipschitz stable for 
a low frequency component of c~ 2 under the assumption that the 
Ricmannian manifold (M 1 c~ 2 dx 2 ) has a strictly convex function 
with no critical points. That is, we show that for all bounded 
C 2 neighborhoods U of c, there is a C 1 neighborhood V of c and 
constants C, R > such that 



\F(c- 2 -c~ 2 ) (0\<Ce 2R ^ 



A- A 



for all c 6 U ("1 V, where A is the Dirichlet-to- Neumann map corre- 
sponding to the wave speed c and 1 1 - 1 1 ^ is a norm capturing certain 
regularity properties of the Dirichlet-to-Neumann maps. 



1. Introduction 

Let M C W 1 be a compact set with nonempty interior and a smooth 
boundary dM and let c G C°°(M) be strictly positive. We consider 
the wave equation on M, 

(1) d?u(t,x) - c(x) 2 Au(t,x) = 0, (t, x) E (0, oo) x M, 

u(t, x) = fit, x), (t, x) e (0, oo) x <9M, 

u|t=o(a?) = 0, d t u\ t =o{x) =0, x G M. 

Let us denote the solution of ([I]) by u*{t,x) = u(t,x) and let T > 0. 
We define the operator 

A CjT : / H- cW|(o,T)x9M, / e C7 CO ((0,T) x 9M), 
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where d u is the Euclidean normal derivative on dM. Often we write 
At '■= A c t- The operator models acoustic boundary measurements 
and is called the Dirichlet-to-Neumann operator. Let us denote by 
C^°(M) the set of the strictly positive functions in C°°(M). Then the 
inverse problem for the wave equation can be formulated as follows: 

(IP) Reconstruct c G C^°(M) given the operator A c> t- 

The finite speed of propagation for the wave equation ([!]) gives a 
necessary condition for T in order to (IP) to have a unique solution. 
Indeed, if there is x G M such that T < 2d(x ,dM), where d is 
the distance function of the Riemannian manifold (M,c~ 2 dx 2 ), then 
the measurements can not contain any information about c{xq). 
Conversely, the problem (IP) is known to be uniquely solvable for T 
strictly greater than the maximum of 2d(x, dM) for x G M. The global 
uniqueness can be proven either by using the Boundary Control (BC) 
method originated from [7] or by using the complex geometric optics 
(CGO) solutions originated from [3D]. However, a typical application 
of the BC method depends on Tataru's unique continuation theorem 
[31], whence only logarithmic type stability is expected for such an ap- 
plication. The CGO solutions based approach is also typically limited 
to logarithmic type stability |3Uj . 

Here we will introduce a global reconstruction method and prove that 
it is locally Lipschitz stable in the sense that we will describe below. 
The method is a modification of the BC method and employs also the 
harmonic exponential functions of the form 

e ,(M - l/2 , ^6 1", |f | = \r)\, £1-77, 

that are CGO solutions for the Euclidean Laplacian. 

Holder stability with an exponent strictly better than 1/2 allows 
an inverse problem to be solved locally by the nonlinear Landweber 
iteration |15| . Moreover, the convergence rate of the iteration is linear 
if and only if the problem is Lipschitz stable. Hence Lipschitz stability 
for (IP) would be important even without our explicit reconstruction 
method. 

Holder type stability results for (IP) were first obtained in [SUEZ], 
and the best Holder exponent available in the literature is 1/2, see 
[9j. However, the Holder exponent 1/2 does not allow the convergence 
result [15] to be applied in a straightforward manner. Moreover, the 
technique in [9J does not give a global reconstruction method since it 
employs the geometric optics solutions corresponding to a fixed wave 
speed c , and requires c to be known a priori. 
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The stability result in [9] depends on the assumption that the Rie- 
mannian manifold (M, CQ 2 dx 2 ) is simple. Similarly, our result depends 
on an assumption that we call stable observability (see Definition 1 be- 
low) and that is also of geometrical nature. Let us also point out that 
Stefanov and Uhlmann [38] have considered the linear inverse problem 
to recover the initial value of a solution u to the wave equation in W 1 
given its trace on dQ x (0,T) where Q C M. n . They prove that recon- 
struction of u(0) supported in Q can not be Lipschitz stable if there 
is a geodesic 7 such that 7(0) G f2 and 7Q— T, T]) n dQ = 0. Thus 
it would be unexpected if the non-linear inverse problem (IP) was lo- 
cally Lipschitz stable without additional assumptions on the geometry 
(M, c- 2 dx 2 ). 

1.1. Statement of the main results. We recall that the wave equa- 
tion ([T]) is said to be continuously observable from open V C dM in 
time T > if there is C f, s > such that 

(2) \\u(T)\\ Hlo(M) + \\d t u(T)\\ L2(M) — Cobs ^\d v U^ L 2^ QT ^ r ^ , 
where u is a solution of the wave equation 

(3) d 2 u - c 2 Au = 0, in (0, T) x M, 

u|(o,t)x9m = 0, in(0,T)x<9M. 

The condition by Bardos, Lebeau and Rauch gives a geometric char- 
acterization of the continuous observability [6j [T3J. In particular, if 
r = dM and (M, c~ 2 dx 2 ) is non-trapping then the continuous observ- 
ability (j2J) is valid. This is analogous with the condition in the above 
mentioned [38] . We refer to [6] for the precise formulation of the geo- 
metric condition. 

Definition 1. Let U C C^°(M). We say that the wave equations (QJ) 
are stably observable for c G U , from open V C dM in time T > 0, if 
there is C b s > satisfying the following: for all c G U the solutions of 
the wave equation (T3|) satisfy the observability inequality (TJ|). 

This stronger form of observability does not follow from the technique 
in [6] since the compactness-uniqueness argument there does not bound 
the constant C & s in terms of the geometry (M,c~ 2 dx 2 ). However, we 
will prove the following theorem. 

Theorem 1. Let c G C^°(M) and suppose that there is a strictly convex 
function i G C 3 (M) with respect to the metric tensor c~ 2 dx 2 , and that 
I has no critical points. Let U C C^°(M) be bounded in C 2 (M) and let 
T C dM be a neighborhood of 

{x G dM; W(x) ■ v > 0}. 
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Then there is a neighborhood V of c in C 1 (M) and T > such that the 
wave equations (QJ) are stably observable for the wave speeds in the set 
U H V, from T in time T. 

We will prove Theorem[T]in Section [4] also for anisotropic wave speeds 
by first deriving a geometric Carleman estimate for the wave equa- 
tion. The main feature of the estimate is the absence of lower order 
terms, and we will follow the tradition of this type of estimates, see 
[2SJ [27J H2J S3] , where the continuous observability was studied but the 
dependence of the constant C b s on the coefficients of the equation was 
not considered; and [16], where the dependence of C Q b s on the zeroth 
order coefficient of the equation was studied. 

There are non-simple Riemannian manifolds (M, g) that admit strictly 
convex functions with no critical points (a trivial example being a 
non-convex subset of the Euclidean space). See [19] and [31] for fur- 
ther discussion on the relations between simplicity, the existence of a 
strictly convex function and the characterization by Bardos, Lebeau 
and Rauch. 

To formulate our main result, let us recall that the Dirichlet-to- 
Neumann operator is continuous, 



A r : F^((0,T) x dM) L 2 ((0,T) x dM), 
where H^((0,T) x dM) := {/ G ^((O.T) x dM); f(0,x) = 0}, 



see [2S] • Moreover, in Section [2J we show that A 2 t has the additional 
regularity-symmetry property, 



where K(A 2 t) '■= RA T RJQ — JA 2 t@, R is the time reversal on (0, T), 
that is Rf(t) := f(T — t), is the extension by zero from (0,T) to 
(0, 2T) and 



The additional regularity can be understood by noticing that K(A) 
corresponds formally to the operator (A* — A) J+ [A, J], where A = A 2 t 
and J are operators of orders -1 and zero, respectively. We denote 
T := (0,T) x dM and define 



(4) 



K(A 2T ) : L 2 ((0,T) x dM) -> L 2 ((0,T) x dM), 




A 



2T 



||-^(A2t)|| L 2 (x) ^ L 2 (t) 



+ I|A< 



T|lHi c (T)^L 2 (T) • 



We will show that 1 1 - 1 1 ^ is a norm in the appendix below. Our main 
result is the following. 
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Theorem 2. Suppose that the wave equations (OP are stably observable 
for the wave speeds in a set U C C+(M), from dM in time T > 0. 
Suppose, furthermore, that there is €u > such that 

c(x) > eu, for all x G M and c EU . 

Let R > satisfy M C B(0,R) and let c e U. Then there is C > 
depending on M, T, c, eu and C b s such that for allcE U 

\F (c~ 2 - c- 2 ) (01 < Ce 2R ^ ||A,, 2T - A c>2T \l , £ G W\ 

where J^ip), p G C°°(M), denotes the Fourier transform of the exten- 
sion by zero of p onto ~R n . 

1.2. Previous literature. We refer to the monograph [23] and to the 

review article [8] for literature concerning the BC method, and to the 
review article [H] concerning the CGO solutions. A stability result 
without a modulus of continuity for the former approach was proved 
in [1] and the first logarithmic type stability result for the latter in [I]. 
Holder type stability results for (IP) are proved in the above mentioned 
articles J^JEZIE]. 

There has been recent interest in results showing that the ill-posedness 
of the inverse problem for the Helmholtz equation decreases when the 
frequency increases, see [22j [32] and the references therein. Moreover, 
in a recent preprint [3] , Lipschitz stability for determining the low fre- 
quency component of a potential in the inverse scattering problem was 
established. This result is similar in spirit with the present one but it 
is based on different techniques. 

As for recovering the potential in a wave equation, Holder type sta- 
bility result was first established in [39J . This was improved later to an 
almost Lipschitz stability result with the Holder exponent being 1 — e in 
[5] . Lipschitz type stability can be obtained if the potential is assumed 
to be parametrized in a finite dimensional space [33] . See also [21 [33] 
for a Lipschitz stability result with finite number of parameters. 

The inverse problem in the present paper is formulated with many 
boundary measurements. On the other hand, for a different formula- 
tion of the inverse problem with a single measurement, Lipschitz type 
stability results can be achieved. For example, [3S] proved Lipschitz 
stability for recovering the sound speed from a single measurement in 
the context of multi-wave imaging. However, such formulation typi- 
cally requires non-vanishing initial data which is not in favor of prac- 
tical applications if only acoustic waves are used for imaging. The 
main methodology used in the inverse problems for the wave equation 
with a single measurement is based on Carleman type estimates. The 
technique was originated in [13] and has been developed tremendously 
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since then. In particular, the continuous observability inequality (EJ) 
may be used to derive the Lipschitz type stability. For more details 
about the single measurement formulation, we refer to [20| l21j [29] and 
the references therein. 

Let us also point out that the continuous observability inequality ([2]) 
is equivalent to the exact controllability of the wave equation, i.e. the 
surjectivity of the control to solution map. This well-known link to the 
control theory has been well studied since 1980s and we refer to the 
review articles [2"8"| ITS] and the references therein for more details on 
the exact controllability of wave equations. 

2. A MODIFICATION OF THE BOUNDARY CONTROL METHOD 

Let p G C°°(M) and let us extend p by zero to M. n . We denote by 
•F(p) the Fourier transform of the extension. Moreover, let us define 
the operator 

B(A CiT ) := RA CtT RI% - IT U 

where Tj, j = 0, 1, are the first two traces on dM, that is %<j) = 4>\dM 
and 7i0 = d u cj)\dM, and 

If(t) := ^ f(s)ds, feL 2 (0,T),te (0, T). 

In this section we will prove the following reconstruction formula. 

Theorem 3. Let c e C™{M), open T C dM and T > satisfy the 
continuous observability inequality (TJ|). Let £, rj G lR n satisfy |£| = \rj\ 
and £ _L rj, and define the functions 

(j)^(x) := e^ +iri) - x/2 , ^ v {x) := e^'^ 2 . 

Then 

(5) 7 (c~ 2 ) (0 = (K(A c , 2T ) t J B(A c , T )^, r) , B{A C , T )^ ,ri)L 2 ((0,T)xdM), 

where K(A c> 2t) is the pseudoinverse of K(A c 2 t) ■ 

We refer to [17] and [10] for the definition and general theory of 
pseudoinverse operators. By inspecting the proof of Theorem [3] we see 
that a formula of the type (jSJ) can be obtained given only the restriction 

A c ,T,r : / H- ^Ifo.Tjxr, / e C™((0,T) x dM), 

of the Dirichlet-to-Neumann map. However, we will leave the partial 
data problems as an object of future study. In particular, the local 
data problem, with the sources / supported on (0, T) x T, will require 
further analysis. 
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2.1. Blagovescenskii type identities. The BC method is based on 
the following identity that originates from |12j . 

(6) (u f (T),u h (T)) L 2 {M . c ~2 {x)dx) = (/, K(A Ct 2T)h) L 2(( 0tT ) xdM ), 

where /, h G Cg°((0,T) x <9M), T > and u f denotes the solution of 
(CQ). The formulation of the identity fl6]) by using the operator K(A c ^t) 
is from [11] . Let us also mention that the iterative time-reversal control 
method introduced there can be adapted to give an efficient implemen- 
tation of the reconstruction formula ([5]). 

We define the map W Ct xf '■= u^(T) and have, sec [26J, 

W CyT : L 2 ((0,T) x DM) -»■ L 2 (M). 

To simplify the notation, we often write K := K(A Ci 2t), B := B(A c t) 
and := W C>T - Then © yields K = W*W, and if extends as a 
continuos operator on L 2 ((0,T) x <9M) by the L? continuity of W. 

Lemma 1. Let f G C^{{0, T) x DM) and let <p G C°°(M) be harmonic. 
Then 

(7) (u f (T), 4>) L 2 {M;c^{x)dx) = (f, B(j)) L 2(( , T )xdM)- 

In particular, B is the restriction of W* on harmonic functions. 
Proof. Let t G (0,T). Then integration by parts gives 

dl{u f {t), (f))L 2 {M;c- 2 {x)dx) 

= (Au f (t),<l)) L 2 {M . dx) - (u f (t),A(j)) L 2 {M . dx) 
= (A T f(t),(j)) L 2( dM) - (f(t),d u (f)) L 2 {dM) . 

By solving this differential equation with vanishing initial conditions 
at t = we get 

(u f (T), (j))L2(M;c- 2 (x)dx) 



(A T f(t),(j)) L 2 {dM) - (f(t),d u (f)) L 2 idM) dtds 

L 2 ((0,T)xdM) — (%f, d u 4>)L 2 ((0,T)xdM), 

where X/(s) := J* f(t)dt. The equation (j7j) follows since = RA T R 
and 1* = I. □ 

2.2. Computation of boundary controls. Let (ft G L 2 (M) and con- 
sider the control equation 

(8) W/ = 0, for / G L 2 ((0, T) x <9M). 

Typically W is not injective and we can hope to solve © only in 
the sense of the pseudoinverse / = W*(j). It is well-known that the 
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pseudoinverse W' is a bounded operator if and only if W has closed 
range R(W). If T > is large enough then R(W) is dense in L 2 (M) 
by Tataru's unique continuation [JT] . Hence the pseudoinverse W ' is a 
bounded operator if and only if W is surjective. It is well-known, see 
e.g. [6], that the map 

/ i — y (u f (T),d t u(T)) : L 2 ((0,T) x dM) -> L 2 (M) x R-\M) 

is surjective if and only if the continuous observability (j2J) holds with 
T = dM. Let us now assume that (J2J) holds. Then 

H/ f : L 2 (M) -> L 2 ((0,T) x 9M) 

is continuous and WW^ is the identity operator. 

The pseudoinverse can be written as = (W*WyW*. Moreover, 
R{K) = R(W*W) = R(W*) is closed since R(W) is closed. Hence 
is continuous on £ 2 ((0, T) x dM) and KK^ is the orthogonal projection 
onto R(K). In particular, if G C°°(M) is harmonic, then W'(j) = 
K^Bip can be computed from the boundary data A 2 t- 

Let <fr,ip E C°°(M) be harmonic. Then the identity (J7|) yields 

(0, 1p)L2(M;c-2(x)dx) = {WW ] (j), i)) L 2(M;c-2(x)dx) 

= (W j (j),B4>) L 2 mT)xdM) = (K ] B(p, Bi)) L 2 mT)xdM) . 

Notice that the functions <fi := <j)^ tV and i/) := ip^ v defined in Theorem 
[3] are harmonic and (f>(x)ip{x) = e l ^' x . Thus 

F { C ~ 2 ) (0 = i. ( P^)L^M ] c^{x)dx) = (K^Bcj), -B'0) i 2(( OiT ) xaM ). 

This proves Theorem [3J 

3. Stability of the reconstruction 

Let us assume that c G C^°(M), T > and T = dM satisfy ©. We 
denote T := (0,T) x dM. Notice that W*<p = d v u\y, where u is the 
solution of the wave equation ([3]) with (u, d t u) = (0, 0) as the initial 
data at t = T. Hence we have 

U\\ L 2 {M ) < Cobs ||W>|| i2(T) . 

In particular, W* is an injection and (W*)~ l : R(W*) -> L 2 (M) satis- 
fies || (W*)- 1 ]]^^^^ < C obs . Notice that 

where Prsw*) is the orthogonal projection onto R(W*). Hence 

11^ lll,2(Af)->.L a (T) = IK^ ) IL 2 (T)-hL 2 (Af) - Cobs- 
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Moreover, (WW)" 1 " = Wt(W*)t, which implies 

11-^ IL 2 (T)->L2(T) = II ^) IL 2 (T)-kL2( T ) — ^obs- 

We are now ready to prove Theorem [2] formulated in the introduction. 

Proof of Theorem® Let c, c G U and denote At = A C) t, and Ay = Ac,t 
and let us define W, W, K, K, B and B analogously. From now on 
we will omit writing L 2 (T) as a subscript. We have, see e.g. 



K ] -K ] 



< 3max( 



A'* 



\mf 



K-K 



K-K 



Notice that for (ft G L7°°(M) the function 

(RI%<}>){t,x) = t<f)(x), te[0,T], xe dM, 

is in H^ C (T). Thus there is Co > depending only on T and M such 
that 



{B-B)cft = {A T -A T )RI%<p 



<C 



A T - A 



T 



Hi c (T)^L2(T) 

Moreover, if G C°°(M) is harmonic then 

ww = \\w*<t>\\ < a 

w 



II^IIC! 



L 2 (M;c- 2 dx) ' 



where we have denoted C c :— \\w \\L 2 (M)^L 2 ny 
Let (ft, ift G C°°(M) be harmonic. Then 

|(0,^)L2(M;?-2dx) - (<P,1p)L 2 (M;c- 2 dx)\ < \((K f ~ K ] )B(ft, Blft) \ 

+ \(K^B(f>, (B-B)ift)\ + \(K^(B-B)(f>,Bi;)\. 

We have 

\((K? -K ] )B(ft, Bift)\ 



K — K 



KK^Bcft, (B-B)tP)\ 



Ar - A. 



T 



L 2 (M;c~ 2 dx) \\Y\\L 2 (M;c- 2 dx) , 



ffi c (THL2(T) 



L 2 (M;c- 2 dx) \\Y\\C 1 (dM) ■ 



Notice that (A^)* = AT f since K* = K. Moreover, as K^Bijj = W^(ft, 
we have 

\{K\B - B)(ft, Bift)\ = \((B- B)<t>,WH)\ 



<C n C, 



O^obs 



A T - A T 



H} C (T)^L 2 (T) 



C x (dM) \\r\\L 2 {M;Z- 2 dx) ■ 



10 



SHITAO LIU AND LAURI OKSANEN 



Hence there is a constant C > depending on M, T, c, eu and C Q b s 
such that for all c £ t/ and harmonic E C 1 (M) 

(9) 

1(0, ^)L 2 (Af;c- 2 da:) ~ (0, ^) ^(M-^dx) \ 



< c 



K - K 



A T - A 



T 



Hi c (T)^L2(T) 



C!(M) 



C X {M) 



Let £ M n satisfy |f| = |r/| and define <j>(x) := e i(?+ir ? W2 . Then 
|0(x)| < e R l«l/ 2 and 

|^-0(x)| < (101 + |^|)e^'/ 2 < 2|e|e^l/ 2 < C R e R ®, 

where Cr > is a constant depending only on R. Hence implies 
that there is C > such that 



|^(c" 2 - c - 2 ) (£)|<Ce 2R KI 



A 



2T 



A 



2T 



£ £ M n , c£ £/. 



□ 



4. Stable observability 

In this section we will prove Theorem [1] formulated in the introduc- 
tion. As the proof is of geometric nature, we will consider the wave 
equation 

d 2 u - A g ^u = 0, in (0, T) x M, 

w|(o,t)x9m = 0, m(0,T)x3M. 

on a smooth compact Riemannian manifold (M, g) with boundary. 
Here A 9iM is the weighted Laplace-Beltrami operator, 

A gifl u = {T l div g (iiVgu), 

where /i £ C°°(M) is strictly positive and div 9 and V g denote the 
divergence and the gradient with respect to the metric tensor g. To 
prove Theorem [1] we will apply the results proven in this section to 
g = c(x)~ 2 dx 2 and n(x) = c(x) n ~ 2 . Then A 9j/i = c(x) 2 A, where A is 
the Euclidean Laplacian. 

We denote by | • \ g , (-, -) g , dV g , dS g , v g and D 2 the norm, the inner 
product, the volume and the surface measures, the exterior unit normal 
vector and the Hessian with respect to g, and will omit writing the 
subscript g when considering a fixed Riemannian metric tensor. We 
recall that the Hessian satisfies 



D 2 w(X, Y) = (D x Vw, Y), D 2 w(X, Y) = D 2 w(Y, X), 
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where w G C 2 (M), X, Y G TM and Dx is the covariant derivative 
with respect to g. We denote div^X := div(/^X) and have the 
formula 

div^wX) = (Vw,X) + wdiv ll X. 

We will obtain a stable observability inequality from a Carleman esti- 
mate similar to that in 1431. 



Lemma 2. Let £ G C 2 (M), $ G C\M) and let w G C 2 (R x M). We 
denote <fi := A^l — if; and q := <fi — |W| 2 and define 

v:=e l w, $:=ipv + 2(Vv } V£), Y := ((d t v) 2 + \Vv\ 2 - qv 2 )V£. 
Then 

(10) e 2£ (<9> - A^wf/2 - d t ($d t v) + div^Vv) + div„ F 

= (d 2 v- A^v + qv) 2 /2 + $ 2 /2 

+ 0(<V) 2 - (j)\Vv\ 2 + 2D 2 £(Vv, Vv) 
+ u(V^, Vv) + (# - div ft (gV^))w 2 . 

Moreover, 

# - div M (gW) = 0| W| 2 + 2Z) 2 £(W, W) 
+ <^-div M (0Vf). 

Proof. Notice that 

A^v = div At (V(e^iu)) = div M (t>W + eVw) 

= (Vv, W) + vA/ + (Ve f , Vic) + e e A„w, 
= 2(Vv, VI) + {A/ - \V£\ 2 )v + e'ApW, 

since 

(Ve^Vw) = (eVl,Vw) = (V£,V(e*w)) - (V£,wVe £ ) 

= (w, Vv) -v\ve\ 2 . 

We expand the square 

(11) e 2l {d 2 w - A^w) 2 = (d 2 v - A^v + qv + $) 2 

= (d 2 v - A^v + qv) 2 + d 2 + 2(d 2 t v - A^v + qv)#, 
and study the cross terms. We have 

(12) W 2 v = dt(&dtv) - dt-d d t v 

= dt(0dtv) - tp(d t v) 2 - 2d t v(Vd t v, V£) 

= d t (W t v) - divJ(d t v) 2 V£) + (A J - i>){d t v)\ 
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since 

2d t v(Vd t v,V£) = (V(d t v) 2 ,V£) = dw^(d t v) 2 V£) - {d t v) 2 A/. 
Similarly, 
(13) 

$A^v = div^tiVv) - (W, Vv) 

= div^VVv) - v(Vip, Vv) - ip\ Vv\ 2 - 2(V(Vt>, VI), Vv), 
= diw^Vv - |Vv| 2 W) - Vt») + (A^-^)|Vt»| 2 

-2D 2 £(Vv, Vv), 

since 

2(V(Vv, W), Vu) = 2{D Vv Vv, VI) + 2(Vu, £> V *V^) 
= 2D 2 v(Vv, V£) + 2D 2 £(Vv, Vv), 

and 

2D 2 v(Vv,V£) = 2D 2 v(V£, Vv) = 2{D m Vv,Vv) = {VI, V\Vv\ 2 ) 
= dw^\Vv\ 2 V£) - \Vv\ 2 A^£. 

Finally, 

(14) §qv = q^v 2 + 2qv(Vv, Vi) 

= (# - div,(qV£))v 2 + div> 2 gW), 

since 

2qv(Vv,V£) = (Vv 2 ,qV£) = div^qVi) - v 2 div ll (qV£). 

The first claim follows by inserting (lT2"j) . ( II 3p and (I14p into (II ip . 
For the second claim notice that 

# - div„(gV*) = (0 - \V£\ 2 )iP - div M ((0 - |W| 2 )W) 
= (A/ - V) | V£\ 2 + 2D 2 £{V£, VI) 
+ #-div M (0VI), 

since 

div M (|W| 2 Vf) = (V|W| 2 , VI) + |W| 2 A/, 

and 

(V|W| 2 , VI) = 2{D vt V£,V£) = 2D 2 £(V£, V£). 

□ 
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Corollary 1 (Pointwise Carleman inequality). Let I G C 3 (M) and 
p > satisfy 

(15) D 2 i(X,X) > p\X\ 2 , X G T X M, x <E M. 
Letr>0 and w G C 2 (R x M). We define 

v := e T V $ := t((A/ - p)v + 2(Vv, W)), 
F := r((c^;) 2 + |Vv| 2 - (rp - r 2 |W|> 2 )W. 

T/ien 

e 2T£ (<9> - A^) 2 / 2 - d t {$d t v) + div M (i9Vw) + div^ Y 

> e 2ri {pr - l){(d t w) 2 + |Vuf )/2 + e 2ri {2p\ W| 2 r - Ci)rV, 
w/iere Ci = p 2 + max x . eA/ |V(A M £(a;))| 2 . 

Proof. We invoke Lemma [2] with £ replaced by t£ and ^ = r(A M £ — p). 
Then <fi = rp. Notice that the two first terms on the right-hand side of 
ffTUj) are positive. We employ (JT5j) for X = Vi> and for X = W to get 

e 2ri (d 2 w - A fl w) 2 /2 - d t {$d t v) + div^Vv) + div M Y 

> pr(d t v) 2 - pr\ Vv\ 2 + 2pr\Wv\ 2 + u(rV(A^ - p), Vv) 

+ (p| W| 2 + 2p| W| 2 )rV + (p(A/ - p) - p div M (W))rV 

> pr(<9^) 2 + (pr - l)|Vw| 2 + 3p|W| W - (p 2 + |VA/| 2 )rV. 
The claim follows by noticing that 

e~ 2re \Vv\ 2 = \twV£ + Vuf > -| Vuf - |W| W. 

□ 

Lemma 3. Let T > and Ze* u; G C 2 ([0,T] x M) safe/y w(t,x) = 
/or (t, x) G [0, T] x dM. Let r, p > and tet £ G C 3 (M). We de/me u, 
$ and Y as in Corollary^ Moreover, we denote dm := fidV and 

(16) r := {x G <9M; (W, i/) > 0}. 
Then 

[ [ div^tiVv) +div l ,Ydmdt < C 2 e BeT r [ I ' (d u w) 2 fidSdt, 

JO J M Jo Jt 

where = 2max xe M^{x) and C 2 = 3vaax xe M \ V£(x)\. Moreover, 

[ \$d t v\dm < (C 3 t + C 4 )e BeT r [ (d t w) 2 + | Vw\ 2 dm, 
Jm Jm 



where 



C 3 = CVmax |W(x)| 2 , C 4 = max |W| + C F max | A J - p|/2 

x£M x£M x£M 
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and Cf > 1 is a constant satisfying the Friedrichs' inequality 
(17) f (j) 2 dm <C F [ |V0| 2 dm, G C °°(M). 

Proof. Notice that on [0, T] x <9M we have v = dtv = and 

= e Tl V«) + rvVi = e Tl Vw = e Tl (Vw, v)v, 
since w vanishes there. Thus 

0(Vv, i/) + (y, i/) = 2r(e r£ (Vu;, V£)(e r£ Vw, v) 
+ r|e r£ (Vw;,z/)z/| 2 (W, i/) 
= 3e 2r£ r(Vu>,i/) 2 (W, v). 

By the divergence theorem, 

/ div(//0Vv + fiY)dV = [ p($Vv + Y,v)dS, 

JM JdM 

and the first claim follows. 

For the second claimed inequality, notice that 

T- 1 e- 7frt \&d t v\ = |((A/ - p)w + 2(Vw + rwVe, V£))d t w\ 

< \{A/- p) +2T\W£\ 2 \\wd t w\+2\W£\\Ww\\d t w\ 

< \(A/ - p)/2 + r\V£\ 2 \(w 2 + (d t w) 2 ) + \V£\(\Vw\ 2 + (d t w) 2 ). 
Hence 



' T r- V I \§d t v\dm 

M 



< (max 



|A/-p|/2 + rmax|W| 2 ) ( / w 2 dm+ / (^w) 2 ^ 
+ max|W| / |Vw| 2 + (d t w) 2 dm, 

JM 

and the second claimed inequality follows from (fTTj) with <fi = w. □ 
Remark 1. Lei w G C 2 ([0, oo) x M) be a solution of 

(18) d 2 u(t,x) — A^u(t,x) = 0, (t, x) G (0, oo) x M, 

u(t,x)=0, (t, x) G (0, oo) x dM, 

Then the energy, 

E(t) := [ (d t u(t)) 2 + |Vu(t)| 2 dm, 



M 



is constant for t G [0, oo). 
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We recall that the constant C\ is defined in Corollary [T] and the 
constants C2, C3, C4 and are defined in Lemma [3J Moreover, we 
define the constant = 2mm xGM £(x). 

Theorem 4 (Observability inequality). Suppose that there is a strictly 
convex function i G C 3 (M) with no critical points. Let p, r > satisfy 

D 2 £(X,X)> p\X\ 2 , \V£(x)\>r, 

for all X G T X M and x G M , and let T C dM contain the set l[Tb)). 
Suppose that 

(19) T > 2(C 3 r + L7 4 )e ( ^- A)r r, w/iere r = max ( -, J . 

\p 2pr l J 

Let u G C 2 ([0,T] x M) be a solution of (El). T/jen 



£(0) < C(T) / / (d v uyndSdt, 
Jo Jv 

where 

Coe {Bi ~ Mr T 

C(T) 



T - 2{C 3 t + CJe&'-MTr ' 
Proof. We will integrate the inequality of Corollary [TJ Notice that 

pr-l>2 and 2p\V£\ 2 r - d > 0. 
By Lemma [3] and Remark [1] 

e ftr / / {d t uf + |Vw| 2 dmdt 

< 2(c7 3 r + C 4 )e B ' T TE(0) + C 2 e BeT r J (d u u) 2 pidSdt. 
To conclude notice that 



[ [ (d t u) 2 + \Vu\ 2 dmdt = TE(0). 
Jo Jm 



□ 



Corollary 2 (Stable observability). Suppose that there is a strictly 
convex function £ G C 3 (M) with no critical points. Let p, r > satisfy 

D 2 £(X,X) > p\X\ 2 , \V£(x)\ > r, 

for all X G T X M and x G M. Suppose that open F C dM satisfies 

{x G dM; (W, v) > 0} C T. 
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Let Uq be a bounded C 2 neighborhood of (g,fj). Then there is a C 1 
neighborhood U of g and constants C, T > satisfying the following: 
for all (g, Jl) G Uq such that g G U, the solutions 

(20) u G C([0,T];H 1 (M)) D C 1 ([0,T];L 2 (M)) 

of the wave equation, 

d 2 u - A~ 9jV u = 0, on (0, T) x M, 

u = 0, on (0, T) x dM, 

satisfy the observability inequality 

Wh^m) + ll^(°) \\l*(m) < C \\ d M\v((o,T)xr) ■ 

Proof. Let us choose a finite number of compact coordinate neighbor- 
hoods covering M and let K be one of them. A metric g is given in 
K by a smooth matrix valued function g^. Let us denote by a(g,x) 

the smallest eigenvalue of the matrix d 2 j£ — T^dki — pgij, where are 
the Christoffel symbols corresponding to "g^. Then a : C 1 x K — > R is 
continuous on the compact set Kq x K, where Kq is the C 1 closure of 
the projection of Uq on the metric tensors. In particular, there is a C 1 
neighborhood U of g such that a(g,x) > in U fl A' x K. That is, 

D|£(X, X) > p\X\i X G TK, geun K . 

The function (g, x) t— >• |Vg£(x)|g is continuous on the compact set K x 
A, whence by making U smaller if necessary, we have 

\Vg£(x)\g >r, x G A, g eU(l K . 

Let us suppose for a moment that A intersects the set dM\T. We may 
assume that dM fl A is given by a defining function F and that v = 
VgF/\VgF\. The function (g,x) h-> (VgA(x), Vg£(x))g is continuous 
on the compact set Ao x A, whence by making U smaller if necessary, 
we have that 

(u(x),Vg£(x))g < 0, x G (A n <9M) \ T, g E U D A . 

By taking the intersection with respect to the finite cover chosen in the 
beginning of the proof, we see that there is a C 1 neighborhood U of g 
such that all g G U D A satisfy the assumptions of Theorem H] with 
the fixed £, p, r and T. 

Let us show next that the constants Cj(g,Jl), j = 1,2,3,4, stay 
bounded in Uq. We may first work locally in a compact coordinate 
neighborhood A as above. Let us denote by \o(g, x) and \ n (g,x) the 
smallest and the largest eigenvalue of gij(x). The functions Ao and A n 
are continuous C° x A — > M, ^(x) is positive definite and U x A is 
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compact in C° x K, where the closure is in C°. In particular, we may 
choose C > so that on Uq 

C'^Xlg < \X\g < C\X\g, X^TK. 

Moreover, we may choose C > so that also the functions \g(x)\ and 
are bounded below by C" 1 and above by C on Uq x K. Hence 
there is Cp > such that the Friedrichs' inequality 

/ <p 2 dm <C F I \V(j)\ 2 gdm, e C °°(M), 

holds for all (g,Jj) e Z7 o- Now it straightforward to see that Cj(g,Ji), 
j = 1, 2, 3, 4, are bounded on U as they can be expressed in coordinates 
using the derivatives of and Jl up to the second order. 

The map (u(0), d t u(0)) \-> <9j,u|(o,T)xr is continuous from Hq(M) x 
L 2 (M) to L 2 ((0, T) x T) by [26J. Thus we may approximate the initial 
data by smooth compactly supported functions and get the observabil- 
ity also for solutions in the energy class ( 120]) . □ 

Theorem CD follows from Corollary [2] by choosing g = c(x)~ 2 dx 2 and 
li{x) = c(x) n ~ 2 . 

Appendix: A linear space for Dirichlet-to-Neumann 

OPERATORS 

Let us consider the operator 

A : #«((0,T) x dM) -> L 2 ((0,2T) x dM), 

and define the map fC(A) := RAtRJQ — J A and the restrictions 

Arf ■— (Af)\( ,T)xdM, A R f := (Af)\(r,2T)xdM- 

Notice that 2J6/(t) = f(s)ds, whence RJOf(0) = and /C(A) is 
well-defined on H^((0,T) x 3M). We define 

ll^ll* := II^-(^)IIl2( X )_ >L 2( T ) + II^t||#i c (t)-^l 2 (t) > 

where we have denoted T := (0, T) x dM. Let us next show that 
is a norm on 

Hh N := {A : H^T) L 2 ((0,2T) x dM); \\A\l < oo}. 

Clearly, is homogeneous and subadditive. Suppose that \\A\\^ = 
and let / e #«(T). Then A T = and 

/•2T-t 

= -2K(A)f(t) = 2JAf(t) = / A R f(s)ds, t e (0, T). 
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By differentiating, we see that ArJ (2T '— £) = for almost all t £ (0, T), 
whence = and we have shown that \\-\\^ is a norm. 
Let us now consider the Dirichlet-to-Neumann operator 

A 2T : f4((0,2T) x dM) -»■ L 2 ((0,2T) x <9M). 

Let _E : H^. c (0,T) — >■ H^ C (0,2T) be an extension operator, that is, 
£/1[o,t] = /■ Then 

A 2T oE: Hl c ((0, T) x dM) L 2 ((0, 2T) x <9M). 

Moreover, causality of the wave equation ([I]) yields (A 2 t o E) T = A T 
and 

(h,1C(A 2T oE)f) = (u h (T),u f (T)) L 2( M ;c- 2 dx) = (h,K(A 2T )f), 

for all /, ft, £ C£°((0, T) x <9M). Hence the embedding A 2T ^ A 2T ° E 
of the Dirichlet-to-Neumann operators to -ff^jv does not depend on the 
choice of the extension operator E. 
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